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We review advances in the theoretical and computational studies of disk winds, jets 
and outflows including: the connection between accretion and jets, the launch of jets from 
magnetized disks, the coupled evolution of jets and disks, the interaction of magnetized young 
stellar objects with their surrounding disks and the relevance to outflows, and finally, the link 
between jet formation and gravitational collapse. We also address the predictions that the theory 
makes about jet kinematics, collimation, and rotation, that have recently been confirmed by high 
spatial and spectral resolution observations. Disk winds have a universal character that may 
account for jets and outflows during the formation of massive stars as well as brown dwarfs. 



1. INTRODUCTION 

The close association of jets and outflows with protostel- 
lar accretion disks is one of the hallmarks of the accretion 
picture of low mass star formation. The most energetic out- 
flow phase occurs during gravitational collapse of a molec- 
ular cloud core - the so-called Class phase - when much 
of its envelope is still raining down onto the forming proto- 
stellar disk and the disk accretion rate is high. Later, in the 
T-Tauri star (TTS) stage, when most of the original core has 
been accreted and the young stellar object (YSO) is being 
fed by lower accretion rates through the suiTounding Kep- 
lerian accretion disk, the high-speed jet becomes optically 
visible. When the disk disappears in the weak-lined TTS 
(WTTS) phase, the jet goes with it. 

The most comprehensive theoretical picture that we have 
for these phenomena is that jets are highly collimated, hy- 
dromagnetic disk winds whose torques efficiently extract 
disk angular momentum and gravitational potential energy. 
Jets also sweep up ambient molecular gas and drive large 
scale molecular outflows. A disk wind was first suggested 
as the origin of jets from accretion disks around black holes 
in the seminal paper by Blandford and Payne (1982, BP82), 
and was soon proposed as the mechanism for protostellar 
jets {Pudritz and Norman, 1983, 1986). 

Several major observational breakthroughs have taken 
place in the study of jets and outflows since PPIV (held in 
1998). Direct, high resolution spectro-imaging and adap- 
tive optics methods discovered the rotation of protostellar 



jets (Bacciotti et al., 2003). These observations also re- 
vealed that jets have an onion-like, velocity structure (with 
the highest speeds being closest to the outflow axis). This 
work provides strong support for the idea that jets originate 
as centrifugally driven MHD winds from extended regions 
of their surrounding disks (see chapter by Ray et al). Re- 
cently, outflows and disks have also been discovered around 
massive stars (see chapter by Arce et al.) as well as brown 
dwarfs (e.g., Bourke et al. 2005) implying that the mecha- 
nism is important across the entire stellar mass spectrum. 

Major advances in the theoretical modeling of these sys- 
tems have also occuiTed, due primarily to a variety of MHD 
computational studies. Simulations now resolve the global 
evolution of disks and outflows, track the interaction of 
disks with central magnetized stars, and even follow the 
generation of outflows during gravitational collapse. These 
studies shows that jets and disks are closely coupled and 
that jet dynamics scales to YSOs of all masses. 

Our review examines the theory of the central engine of 
jets and its exploration through the use of computer simu- 
lations. We focus mainly on developments since PPIV and 
refer to the review by Konigl and Pudritz (2000, KPOO) for a 
discussion of the earlier literature, as well as Pudritz (2003) 
and Heyvaerts (2003) for more technical background. We 
first discuss the basic theory of disk winds and their kine- 
matics. We then switch to computational studies of jets 
from accretion disks treated as boundary conditions, as well 
as global simulations including the disk. We then examine 
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the innermost regions of the disk where the stellar magne- 
tosphere interacts with the disk, as well as the surface of 
the star that may drive an accretion-powered outflow. Fi- 
naUy, we discuss how outflows are generated during the 
early stages of the gravitational coUapse. 

2. THEORY OF DISK WEVDS 

An important insight into the nature of the engine for 
jets can be gleaned from the observed ratio of the momen- 
tum transport rate (or thrust) carried by the CO molecular 
outflow to the thrust that can be provided by the bolomet- 
ric luminosity of the central star (e.g., Cabrit and Bertout, 
1992); 

Foutflow/Frad = 250{Lbol/lO^ Lq)-°-\ (1) 

This relation has been confirmed and extended by the analy- 
sis of data from over 390 outflows, ranging over six decades 
up to IO^Lq in stellar luminosity (Wu et al, 2004). It sug- 
gests that jets from both low and high mass systems are 
probably driven by a single, non-radiative, mechanism. 

Jets are observed to have a variety of structures and 
time-dependent behaviour - from internal shocks and mov- 
ing knots to systems of bow shocks that suggest long- 
time episodic outbursts. They show wiggles and often 
have cork-screw like structure, suggesting the presence ei- 
ther of jet precession, the operation of non-axisymmetric 
kink modes, or both. Given the highly nonlinear be- 
haviour of the force balance equation for jets (the so-called 
Grad-Shafranov equation), theoretical work has focused on 
tractable and idealized time-independent, and axisymmetric 
or self-similar models (e.g., BP82) of various kinds. 

2.1. Conservation Laws and Jet Kinematics 

Conservation laws are the gold standard in physics, and 
play a significant role in understanding astrophysical jets. 
This is because whatever the details (e.g. the asymptotics, 
the crossing of critical points, the way that matter is loaded 
onto field lines within the disks, etc.), conservation laws 
strongly constrain the flux of mass, angular momentum, and 
energy. What cannot be constrained by these laws will de- 
pend on the general physics of the disks such as on how 
matter is loaded onto field lines. 

Jet dynamics can be described by the time-dependent, 
equations of ideal MHD. The evolution of a magnetized, 
rotating system that is threaded by a large-scale field B in- 
volves (i) the continuity equation for a conducting gas of 
density p moving at velocity v (which includes turbulence); 
(ii) the equation of motion for the gas which undergoes 
pressure (p), gravitational (with potential <&), and Lorentz 
forces; (iii) the induction equation for the evolution of the 
magnetic field in the moving gas where the current density 
is j = (c/47r)V X B; (iv) the energy equation, where e is 
the internal energy per unit mass; and, (v) the absence of 
magnetic monopoles. These are written as: 

% + V.{pv) = (2) 



p(|^ + (v.V)v^ +Vp + pV$-^^ = (3) 
— -Vx(vxB) = (4) 

p(|^ + (v.V)e^ +p(V.v) = (5) 
V.B = (6) 

We specialize to the restricted case of stationary, as well 
as 2D (axisymmetric) flows, from which the conservation 
laws follow. It is useful to decompose vector quantities 
into poloidal and toroidal components (e.g. magnetic field 
B = Bp +5^60). In axisymmetric conditions, the poloidal 
field Bp can be derived from a single scalar potential a(r, z) 
whose individual values, a = const, define the surfaces of 
constant magnetic flux in the outflow and can be specified 
at the surface of the disk (e.g., Pelletier and Pudritz, 1992; 
PP92). 

Conservation of mass and magnetic flux along a field 
Une can be combined into a single function k that is called 
the "mass load" of the wind which is a constant along a 
magnetic field line; 

pvp = fcBp. (7) 

This function represents the mass load per unit time, per 
unit magnetic flux of the wind. For axisymmetric flows, its 
value is preserved on each ring of fleld lines emanating from 
the accretion disk. Its value on each field line is determined 
by physical conditions - including dissipative processes - 
near the disk surface. It may be more reveaUngly written as 

^^""^ = = 

where dM^ is the mass flow rate through an annulus of 
cross-sectional area dA through the wind and d'i is the 
amount of poloidal magnetic flux threading through this 
same annulus. The mass load profile is a function of the 
footpoint radius ro of the wind on the disk. 

The toroidal field in rotating flows derives from the in- 
duction equation; 

B4, = ^{v^-^or), (9) 

where fio is the angular velocity of the disk at the mid- 
plane. This result shows that the strength of the toroidal 
field in the jet depends on the mass loading as well as the jet 
density. Denser winds should have stronger toroidal fields. 
We note however, that the density does itself depend on the 
value of k. Equation (9) also suggests that at higher mass 
loads, one has lower toroidal field strengths. This can be 
reconciled however, since it can be shown from the conser- 
vation laws (see below) that the value of k is related to the 
density of the outflow at the Alfven point on a field line; 
k = (/9^/47r)^/^ (e.g., PP92). Thus, higher mass loads cor- 
respond to denser winds and when this is substituted into 
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equation (9), we see that this also implies stronger toroidal 
fields. We show later that jet coUimation depends on hoop 
stress through the toroidal field and thus the mass load must 
have a very important effect on jet colUmation (Section 3.2). 

Conservation of angular momentum along each field 
fine leads to the conserved angular momentum per unit 
mass; 

1(a) = rvA, — — r = const. (10) 

The form for I reveals that the total angular momentum is 
carried by both the rotating gas (first term) as well by the 
twisted field (second term), the relative proportion being 
determined by the mass load. 

The value of l{a) that is transported along each field line 
is fixed by the position of the Alfven point in the flow, where 
the poloidal flow speed reaches the Alfven speed for the 
first time (ttia — 1). It is easy to show that the value of the 
specific angular momentum is; 



l{a) = fiori = {rA/rof lo- 



in) 



where Iq = VKfira = ^o^q is the specific angular momen- 
tum of a Keplerian disk. For a field line starting at a point ro 
on the rotor (disk in our case), the Alfven radius is rA(ro) 
and constitutes a lever arm for the flow. The result shows 
that the angular momentum per unit mass that is being ex- 
tracted from the disk by the outflow is a factor of (rA/rg)^ 
greater than it is for gas in the disk. For typical lever arms, 
one particle in the outflow can carry the angular momentum 
of ten of its fellows left behind in the disk. 

Conservation of energy along a field line is expressed 
as a generalized version of Bernoulh's theorem (this may be 
derived by taking the dot product of the equation of motion 
with Bp). Thus, there is a specific energy e{o.) that is a con- 
stant along field lines, which may be found in many papers 
(e.g., BP82 and PP92). Since the terminal speed Vp = Voc 
of the disk wind is much greater than its rotational speed, 
and for cold flows, the pressure may also be ignored, one 
finds the result: 



Woo - = {rA/ro)vc 



(12) 



There are three important consequences for jet kinemat- 
ics here; (i) that the terminal speed exceeds the local escape 
speed from its launch point on the disk by the lever arm ra- 
tio; (ii) the terminal speed scales with the Kepler speed as a 
function of radius, so that the flow will have an onion-Uke 
layering of velocities, the largest inside, and the smallest on 
the larger scales, as seen in the observations; and (iii) that 
the terminal speed depends on the depth of the local grav- 
itational well at the footpoint of the flow - implying that it 
is essentially scalable to flows from disks around YSOs of 
any mass and therefore universal. 

Another useful form of the conservation laws is the com- 
bination of energy and angular momentum conservation to 
produce a new constant along a field line (e.g., PP92); 
j{a) = e{a) — fJo^a). This expression has been used (An- 
derson et al, 2003) to deduce the rotation rate of the launch 



region on the Kepler disk, where the observed jet rotation 
speed is ti^^oo at a radius Too and which is moving in the 
poloidal direction with a jet speed ofvp^^o- Evaluating j for 
a cold jet at infinity and noting that its value (calculated at 
the foot point) is j(ao) = —(3/2)?;^ g, one solves for the 
Kepler rotation at the point on the disk where this flow was 
launched: 



/(2t;,^,cx>r'oo)- 



(13) 



When applied to the observed rotation of the Large Veloc- 
ity Component (LVC) of the jet DG Tau (Bacciotti et al., 
2002), this yields a range of disk radii for the observed ro- 
tating material in the range of disk radii, 0.3^ AU, and the 
magnetic lever arm is VA/ro — 1.8-2.6. 

2.2. Angular Momentum Extraction 

How much angular momentum can such a wind extract 
from the disk? The angular momentum equation for the 
accretion disk undergoing an external magnetic torque may 
be written: 



,V d{roVo) _ Or r 
Ma : = —rnO^tS 



dro 



(14) 



where we have ignored transport by MRl turbulence or spi- 
ral waves. By using the relation between poloidal field and 
outflow on the one hand, as well as the link between the 
toroidal field and rotation of the disk on the other, the angu- 
lar momentum equation for the disk yields one of the most 
profound scaling relations in disk wind theory - namely - 
the link between disk accretion and mass outflow rate (see 
KP00,PP92 for details): 



Ma ~ (rA/ro)2M^ 



(15) 



The observationally well known result that in many sys- 
tems, Mw/ Ma — 0.1 is a consequence of the fact that lever 
arms are often found in numerical and theoretical work to 
be ta/^'o — 3 - the observations of DG Tau being a per- 
fect example. Finally, we note that the angular momentum 
that is observed to be carried by these rotating flows (e.g. 
DG Tau) is a consistent fraction of the excess disk angular 
momentum - from 60-100% (e.g., Bacciotti, 2004), which 
is consistent with the high extraction efficiency discussed 
here. 

2.3. Jet Power and Universality 

These results can be directly connected to the observa- 
tions of momentum and energy transport in the molecu- 
lar outflows. Consider the total mechanical power that is 
carried by the jet, which may be written as (e.g., Pudritz, 
2003); 



-'jet 



\ / dM^vl 
J n 



GM^M, 




This explains the observations of Class outflows wherein 
Lw/Lboi — 1/2, since the main luminosity of the central 
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source at this time is due to accretion and not nuclear re- 
actions. (The factor of 1/2 arises from the dissipation of 
some accretion energy as heat at the iimer boundary). The 
ratio of wind to stellar luminosity decreases at later stages 
because the accretion luminosity becomes relatively small 
compared to the bolometric luminosity of the star as it nears 
the ZAMS. 

This result states that the wind luminosity taps the gravi- 
tational energy release through accretion in the gravitational 
potential of the central object - and is a direct consequence 
of Bernoulli's theorem. This, and the previous results, im- 
ply that jets may be produced in any accreting system. The 
lowest mass outflow that has yet been observed corresponds 
to a proto-brown dwarf of luminosity ~ 0.09Lq, a stel- 
lar mass of only 20 — 45Mj„p, and a very low mass disk 
< 1O-^M0 {Bourke etal, 2005). 

It should be possible, therefore, for lower luminosity jets 
to be launched from the disks around Jovian mass plan- 
ets. Recent hydrodynamical simulations of circumstellar 
accretion disks containing and building up an orbiting pro- 
toplanetary core have numerically proven the existence of 
a circum-planetary sub-disk in almost Keplerian rotation 
close to the planet {Kley et al, 2001). The accretion rate 
of these sub-disks is about M^-p = 6 x 10~^ Mj^p yr~^ 
and is confirmed by many independent simulations. With 
that, the circum-planetary disk temperature may reach val- 
ues up to 2000 K indicating a sufficient degree of ionization 
for matter-field coupUng and would also allow for strong 
equipartition field strength (Fendt, 2003). 

The possibility of a planetary scale MHD outflow, sim- 
ilar to the larger scale YSO disk winds, is indeed quite 
likely because: (i) the numerically established existence of 
circum-planetary disks is a natural feature of the formation 
of massive planets; and (ii) the feasibiUty of a large-scale 
magnetic field in the protoplanetary environment (Quillen 
and Trilling, 1998; Fendt, 2003). One may show, moreover, 
that the outflow velocity is of the order of the escape speed 
for the protoplanet, at about 60 km (Fendt, 2003). 

On very general grounds, disk winds are also hkely to 
be active during massive star formation (e.g., Konigl 1999). 
Such outflows may already start during the early collapse 
phase when the central YSO still has only a fraction of a 
solar mass (e.g., Pudritz and Banerjee, 2005). Such early 
outflows may actually enhance the formation of massive 
stars via disk accretion by punching a hole in the infalling 
envelop and releasing the building radiation pressure (e.g., 
Krumholz et al., 2005). 

2.4. Jet CoUimation 

In the standard picture of hydromagnetic winds, collima- 
tion of an outflow occurs because of the increasing toroidal 
magnetic field in the flow resulting from the inertia of the 
gas. Beyond the Alfven surface, equation (8) shows that 
the ratio of the toroidal field to the poloidal field in the 
jet is of the order B^/Bp ~ r/rx '> 1, so that the field 
becomes highly toroidal. In this situation, colUmation is 



achieved by the tension force associated with the toroidal 
field which leads to a radially inwards directed component 
of the Lorentz force (or " z-pinch"); FLorentz,r — jzB^. The 
stability of such systems is examined in the next section. 

In Heyvaerts and Norman (1989) it was shown that two 
types of solution are possible depending upon the asymp- 
totic behaviour of the total current intensity in the jet; 

7 = 27r /" i,(r', z')dr' = (c/2)rB^. (17) 
Jo 

In the limit that 7 — » as r — > oo, the field Unes are 
paraboloids which fill space. On the other hand, if the cur- 
rent is finite in this limit, then the flow is colhmated to cyUn- 
ders. The colUmation of a jet therefore depends upon its 
current distribution - and hence on the radial distribution of 
its toroidal field - which, as we saw earlier, depends on the 
mass load. Mass loading therefore must play a very impor- 
tant role in controlling jet coUimation. 

It can be shown (Pudritz et al., 2006; PRO) that for a 
power law distribution of the magnetic field in the disk, 
BzirojO) oc Tq"^ and an injection speed at the base of 
a (polytropic) corona that scales as the Kepler speed, that 
the mass load takes the form k oc r~^~>^. In this regime, 
the current takes the form 7(r, 2) oc ro^~^^^^\ Thus, 
the current goes to zero for models with /x < —1/2, and 
that these therefore must be wide angle flows. For models 
with /z > —1/2 however, the current diverges, and the flow 
should coUimate to cyhnders. 

These results predict that jets should show different de- 
grees of coUimation depending on how they are mass loaded 
(PRO). As an example, neither the highly centraUy con- 
centrated, magnetic field Unes associated with the initial 
split-monopole magnetic configuration used in simulations 
by Romanova et al. (1997), nor the similar field structure 
invoked in the X-wind (see review by Shu et al., 2000) 
should become coUimated in this picture. On the other 
hand, less centraUy (radiaUy) concentrated magnetic con- 
figurations such as the potential configuration of Ouyed and 
Pudritz (1997a, OPI) and BP82 should colUmate to cyUn- 
ders. 

This result also explains the range of coUimation that is 
observed for molecular outflows. Models for observed out- 
flows fall into two general categories: the jet-driven bow 
shock picture, and a wind-driven shell picture in which the 
molecular gas is driven by an underlying wide-angle wind 
component such as given by the X-wind (see review by 
Cabrit et al, 1997). A survey of molecular outflows by 
Lee et al. (2000) found that both mechanisms are needed in 
order to explain the fuU set of systems observed. 

Finally, we note that apart from these general theorems 
on colUmation, Spruit et al. ( 1 997) has proposed that a suffi- 
ciently strong poloidal field that is external to the flow could 
also force its coUimation. According to this criterion, such 
a poloidal field could coUimate jets provided that the field 
strength decreases no slower than B^ ~ r~'^ with < 1.3. 
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3. SIMULATIONS: DISKS AS JET ENGINES 

Computational approaches are necessary if we are to 
open up and explore the vast spaces of solutions to the 
highly nonlinear jet problem. The first simulations of non- 
steady MHD jets from accretion disks were performed by 
Uchida and Shibata (1985) and Shibata and Uchida (1986). 
These early simulations were based on initial states that 
were violently out of equilibrium and demonstrated the role 
of magnetic fields in launching and accelerating jets to ve- 
locities of order of the Keplerian velocity of the disk. 

The published simulations of non-radiative ideal MHD 
YSO jets differ in their assumed initial conditions, such as 
the magnetic field distribution on the disk, the conditions in 
the plasma above the disk surfaces, the state of the initial 
disk corona, and the handling of the gravity of the central 
star. Nevertheless, they share common goals: to establish 
and verify the four important stages in jet evolution, namely, 
(i) ejection; (ii) acceleration; (iii) collimation; and (iv) sta- 
bility. In the following, we describe a basic physical ap- 
proach for setting up clean numerical simulations and how 
this leads to advances in our understanding of disk winds 
and outflows. 

3.1. 2-Diinensional Simulations 

The simulations we discuss here have been studied in 
greater detail in Ouyed and Pudritz (1997a; OPI); Ouyed 
and Pudritz (1997b; OPII); Ouyed et al. (1997; OPS); 
Ouyed and Pudritz (1999; OPIII); and PRO. We use these 
for pedagogical purposes in this review. To see anima- 
tions of the simulations presented here as well as those per- 
formed by other authors, the interested reader is directed 
to the "animations" link at ' http://www.capca.ucalgary.ca' 
(hereafter CAPCA). These simulations were run using the 
ZEUS (2-D, 3-D and MP) code which is arguably the best 
documented and utilized MHD code in the literature (Stone 
and Norman, 1992, 1994). It is an explicit, finite difference 
code that runs on a staggered grid. The equations generally 
solved are those of ideal MHD listed in Section 1, equa- 
tions (2)-(6). 

The evolution of the magnetic field (induction equa- 
tion (4) above) is followed by the method of constrained 
transport. In this approach, if V.B = holds for the initial 
magnetic configuration, then it remains so for all later times 
to machine accuracy. The obvious way of securing this con- 
dition is to use an initial vector potential A(r, z, < — 0) that 
describes the desired initial magnetic field at every point in 
the computational domain. 

To ensure a stable initial state that allows for a tractable 
simulation, three simple rules are useful (OPS and OPI): (i) 
use a corona that in hydrostatic balance with the central ob- 
ject (to ensure a perfectly stable hydrostatic equilibrium of 
the corona, the point-mass gravitational potential should be 
relocated to zone centers - see Section 4.2 in OPI for de- 
tails); (ii) put the disk in pressure balance with the corona 
above it, and (iii) use a force-free magnetic field configu- 
ration to thread the disk and corona. The initial magnetic 



configurations should be chosen so that no Lorentz force is 
exerted on the initial (non-rotating) hydrostatic corona de- 
scribed above. 

The BP82 self-similar solution for jets uses a simple 
polytropic equation of state, with an index 7 = 5/3. Their 
solution is possible because this choice eliminates the many 
complications that arise from the the energy equation, while 
preserving the essence of the outflow problem. This Ansatz 
corresponds to situations where the combined effects of 
heating and cooling simulates a tendency toward a locally 
constant entropy. This choice also greatly simplifies the nu- 
merical setup and allows one to test the code against ana- 
lytic solutions. 

Another key simplification in this approach is to exam- 
ine the physics of the outflow for fixed physical conditions 
in the accretion disk. Thus, the accretion disk at the base 
of the corona - and in pressure balance with the overlying 
atmosphere - is given a density profile that it maintains at 
all times in the simulation, since the disk boundary condi- 
tions are applied to the "ghost zones" and are not part of 
the computational domain. In part, this simplification may 
be justified by the fact that typically, accretion disks will 
evolve on longer time scales than their associated jets. 

The hydrostatic state one arrives at has a simple ana- 
lytic solution which was used as the initial state for all of 
our simulations and was adopted and further developed by 
several groups. These studies include e.g. low plasma-/? 
monopole-like field distributions (Romanova et al., 1997); 
a magneto-gravitational switch mechanism (Meier et al., 

1997) ; the disk accretion-ejection process (Kudoh et al., 

1998) ; the interrelation between the grid shape and the de- 
gree of flow collimation (Ustyugova et al, 1998); a self- 
adjusting disk magnetic field inclination (Krasnopolsky et 
al, 1999, 2003); the interrelation between the jet's turbulent 
magnetic diffusivity and collimation (Fendt and Cemeljic, 
2002); a varation of the disk rotation profile (Vitorino et 
al, 2002); or dynamo maintained disk magnetic fields (von 
Rekowski and Brandenburg, 2004). More elaborate setups 
wherein the initial magnetic field configuration originates 
on the surface of a star and connects with the disk have 
also been examined (Section 5 - e.g., Hayashi et al, 1996; 
Goodson et al, 1997; Fendt and Elstner, 1999; Keppens 
and Goedbloed, 2000). 

The lesson from these varied simulations seems to be 
that all roads lead to a disk field. Thus, the twisting of a 
closed magnetic field that initially threads the disk beyond 
the co-rotation radius rapidly inflates the field and then dis- 
connects it from the star, thereby producing an open disk- 
field line. Likewise, simulations including dynamo gener- 
ated fields in the disk lead to a state that resembles our ini- 
tial setup (see Section 4). Thus, from the numerical point 
of view, a "fixed-disk" simulation is general and useful. 

The setup requires five physical quantities to be speci- 
fied at all points of the disk surface at all times (see OPI 
for all details). These are the disk density piro); compo- 
nents of the vertical and toroidal magnetic field, Bz{ro) 
and Bfftiro); and velocity components in the disk, w^(ro) 
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Fig. 1. — Left panels: initial magnetic field configurations for 
winds with fi = (OPI), fi = -0.25 (BP), = -0.5 (PP), and 
H = —0.75 (steep) Right panels: final magnetic field configura- 
tions (at t = 400) for each case, with Alfven points (filled circles) 
and fast magneto-sonic points (stars) marked. Note the more open 
magnetic - and streamline - structures as fj, goes down. [Adapted 
from PRO]. 



and ^^(ro), where rg is the radius (cyhndrical coordinates 
are adopted with the disk located at z = 0; rg = r{z — 0)). 
The remaining field component (^o ) is determined by the 
solenoidal condition, while the radial inflow speed through 
the disk is neglected (vr {ro) ~ 0) since it is far smaller than 
the sound speed in a real disk. The model is described by 
five parameters defined at the inner disk radius rj, three of 
which describe the initial corona. The two additional pa- 
rameters describe the disk physics and are i/i, which scales 
the toroidal field in the disk = h'i x {r-Jro), and the 
(subsonic) injection speed of material from the disk into the 
base of the corona, Winj ~ Vz{ro)/v^{rQ) ~ 0.001. Sim- 
ulations were typically run with (500 x 200) spatial zones 
to resolve a physical region of (SOrj x 20ri) in the z and 
r directions, respectively. A resolution of 10 zones per rj 
provides enough dynamical range to accurately follow the 
smooth acceleration above the disk surface. The simula- 
tions were run up to 400ti (where ti is the Kepler time for 
an orbit at the inner edge of the disk). 

A series of magnetic configurations is shown in Fig. [2 
where the disk field is modeled as Bz{ro, 0) oc t'q"^. The 
initial configurations range from the rather well collimated 
(such as the potential configuration of OP97, n = 0.0; and 
the Blandford- Payne configuration, fj, ~ —0.25), to the ini- 
tially more open configurations of Pelletier-Pudritz (PP92; 
/X = —0.5), and steeper (/i ~ —0.75). 

A good example of a configuration that evolves into a 



stationary jet was studied by OPI and is shown in the up- 
per panel of Fig. The simulation shows the existence 
of an acceleration region very close to the disk surface. 
The acceleration from the disk occurs by a centrifugal ef- 
fect whereby, at some point along sufficiently inclined field 
lines, centrifugal force dominates gravity and gas is flung 
away like beads on a wire. Thus, a toroidal field component 
is created because the field lines corotate with the under- 
lying disk. The inertia of the matter in the flow region ulti- 
mately forces the field to fall behind the rotation of the disk, 
which produces the toroidal field component. More pre- 
cisely, beyond the Alfven surface (shown as dots in the right 
hand panels in Fig.Q, the hoop stress induced by the self- 
generated Bcji eventually dominates, which provides the col- 
limation. The ratio of the toroidal to poloidal magnetic field 
along illustrative field lines (e.g.. Figs. 3 in OPS) clearly 
shows that the predominant magnetic field in jets beyond 
their fast magnetosonic (FM) surfaces is the toroidal field 
component. The gas is eventually collimated into cylinders 
parallel to the disk's axis starting at the Alfven surface to 
distances much beyond the FM surface (see Fig.Q. 

The final velocities achieved by such winds are of the 
order of 2 times the Kepler velocity, along a given field line 
(e.g.. Fig. 3 in OPS), which translates to roughly 100-300 
km s^^ for a standard YSO. The general trend is that the 
jet solutions, dominated mainly by the poloidal kinetic en- 
ergy, are very efficient in magnetically extracting angular 
momentum and energy from the disk, as confirmed by sim- 
ulations performed by other groups. In 1000 years, for ex- 
ample, and for a standard YSO, the simulations imply that 
the disk winds can carry a total energy of 3 x lO"'^ ergs, 
sufficient to produce the observed molecular outflows 

In OPII an initial vertical field configuration, was used. 
The simulations show that the strong toroidal magnetic field 
generated recollimates the flow towards the disk's axis and, 
through MHD shocks, produces knots. The knot generator 
is located at a distance of about z ~ Srj from the surface 
of the disk (OPS). Knots propagate down the length of the 
jet at speeds less than the diffuse component of the outflow. 
The knots are episodic, and are intrinsic to the jet and not 
the accretion disk, in this calculation. 

A different initial state was used by Krasnopolsky et al. 
(1999), who introduce and maintain throughout the simula- 
tion, a strong outflow on the outflow axis. Otherwise, they 
choose an initial disk field distribution that is the same as 

3/2 

PP92, a mass flux density pvz oc To , so that k — const.. 
The values of Br and at the disk surface were not fixed 
in their simulations. (This does not guarantee that the disk 
or the boundary at the base will remain Keplerian over 
time). The fixed B^ f^)^ profile adopted in the OP setup 
ensures exactly that - see the discussion around Eq (3.46) 
in OPI). To evade the problems with the initial setup men- 
tioned above, these authors continuously launch a cylindri- 
cal wind on the axis. This imposed jet introduces currents 
which could significantly affect the stability and collima- 
tion of the disk wind. Their disk mass loading would pre- 
dict that the disk wind should not be well collimated (see 
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Section 2.4), whereas their simulation does appear to colli- 
mate. This suggests that their on-axis jet may be playing a 
significant role in the simulation results. 

3.2. The role of mass loading 

The mass-load k, defined in equation (8), can be estab- 
lished by varying; the coronal density profile while keep- 
ing pvp constant {Anderson et ai, 2005), the disk rota- 
tional profile since Vp = VinjVcf, (e.g., Vitorino et ai, 2002), 
the distribution of the poloidal magnetic field on the disk 
(PRO), or the distribution of both the disk magnetic field 
and the mass flow profile (Fendt, 2005). 

The prediction that the mass load determines the coUi- 
mation of the jet was tested in simulations by PRO. Fig. 1 
shows that simulations with fj, = 0, —0.25 coUimate into 
cylinders, while the fj, = —0.5 case (PP92) transitions to- 
wards the wide angle flow seen in the = —0.75 simula- 
tion. These results confirm the theory laid out in Section 
2.4. We also note that each of the simulations mentioned 
above results in a unique rotation profile of the jet that might 
in principle be observable. Going from the potential case to 
the Pelletier-Pudritz configuration, the radial profiles of the 
rotational velocity of the jets scale as power laws 

t;^(r,oo) ocr" (18) 

where a = —0.76; —0.66; —0.46 respectively. The ob- 
servation of a rotational profile would help pick out a unique 
mass loading in this model (PRO). 

The mass loading also affects the time-dependent be- 
haviour of jets. As the mass load is varied a transition 
from stationary, to periodic and sometimes to a discontin- 
uous dying jet occurs. For example, in some of the sim- 
ulations it was found that low mass loads for jets lead to 
rapid, episodic behaviour while more heavily mass-loaded 
systems tend to achieve stationary outflow configurations 
(see OPlll). 

It is interesting that the simulations of Anderson et al. 
(2005) did not find any non-steady behavior for low mass 
loading (differently defined than equation (8)) but instead 
find an instability for high mass loading. They suggest 
that the origin of the non-steady behavior for high mass 
loading is that the initially dominant toroidal field they im- 
pose could be subject to the kink instability (e.g., Cao and 
Spruit, 1994). The very large mass-loads, and large injec- 
tion speeds (winj = (0.01, 0.1)) they use drive an instability 
probably related to excessive magnetic braking. Some of 
the main differences with OP are: (i) the use of the non- 
equilibrium set up of Krasnopolsky et al. which introduces 
a strong current along the axis that could strongly affect the 
stability properties of these outflows; and (ii) the large in- 
jection speeds make the sonic surface too close to the disk 
(the condition that Cg <C i^k is not satisfied) and does not 
provide enough dynamical range for the gas launched from 
the disk to evolve smoothly. Numerical instabilities remi- 
niscent of what was found by Anderson et al. (2005) were 
observed by OP and these often disappear with high enough 



resolution (i.e. dynamical range). We suspect that this type 
of instability will disappear when a proper disk is included 
in the simulations (e.g., Casse and Keppens, 2002). 

3.3. 3-Dimensional simulations 

The stability of jets is one of the principal remaining 
challenges in the theory. It is well known that the purely 
toroidal field configurations that are used to help confine 
static, 3-D, tokamak plasmas are unstable (e.g., Roberts, 
1967; Bateman, 1980). The resulting kink, or helical 
(m = 1) mode instability derived from a 3-D linear sta- 
bility analysis is powered by the free energy in the toroidal 
field, namely B^/Stt (Eichler, 1993). Why are real 3-D jets 
so stable over great distances in spite of the fact that they 
are probably threaded by strong toroidal fields? 

The 3D simulations needed to investigate the importance 
of kink modes are still rare. Early attempts include Lucek 
and Bell (1996) and Hardee and Rosen (1999; see also ref- 
erences therein) who performed 3-D simulations of "equi- 
librium" jets, and found that these uniform, magnetized 
jet models remain Kelvin-Helmholtz (K-H) stable to low- 
order, surface helical and elliptical modes (m = 1, 2), pro- 
vided that jets are on average sub-Alfvenic. This is in ac- 
cord with the prediction of linear stabiUty analysis. How- 
ever, most configurations for jet simulations use rather ad 
hoc prescriptions for the initial toroidal field configuration 
so that it is difficult to assess how pertinent the results are 
to the case of a jet that establishes its own toroidal field as 
the jet is accelerated from the accretion disk. In general, the 
available analytic and numerical results for the stability of 
simple jets show that the fastest growing modes are of K- 
H type. These K-H instabilities are increasingly stabilized 
for super- Alfvenic jets, as Mfm is increased much beyond 
unity. It is also generally known that sub-Alfvenic jets are 
stable. Taken together, these results suggest that 3-D jets 
are the most prone to K-H instabilities a bit beyond their 
Alfven surface, a region wherein their destabilizing super- 
Alfvenic character cannot yet be offset by the stabiUzing 
effects engendered at large super FM numbers. 

There have, as yet, only been a few attempts to simu- 
late 3D disk winds. Contrary to what may be intuitive, it 
is inadvisable to perform these 3-D simulations in cylindri- 
cal coordinates. For one thing, special treatment must be 
given to the "wedge zones" that abut the z-axis (no longer a 
symmetry axis in 3D), and velocities that pass through the 
2;-axis pose a very difficult numerical problem. Secondly, 
even with such technical details in hand, plane waves are 
badly disrupted on passing through the z-axis, and this pro- 
vides an undesirable bias to what should be an unbiased (no 
axis should be preferred over another) 3D calculation. For- 
tunately, using Cartesian coordinates to simulate cylinders 
is feasible with careful setups as has been demonstrated in 
Ouyed et al. (2003; OCP); and Ouyed (2003). These au- 
thors used a 3D version of the OPI setup. 

The central result of this study is that jets survive the 
threatening non-axisymmetric kink (m = 1) mode. The 
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Fig. 2. — Cylindrical Jets - upper panels show snapshots of 20 
magnetic field lines of the 3D simulations performed in OCP, ** 
att — 50, 130, 210, and 240. The two central magnetic field lines 
(dotted lines) originate on the central compact object (illustrated 
by the semi-sphere to the left). The disk axis is along the diagonal 
of the frame (on a 45° angle). Notice the violent jet behavior 
in the in-between frames when the kink mode appears. The jet 
eventually aligns itself with the original disk's axis acquiring a 
cylindrically collimated shape. Corkscrew Jets - The lower panel 
is a map of the jet column density at t — 320 for a simulation 
with a profile, imposed at the accretion disk, that is different 
from the simulation shown in the upper four panels (see OCP). 
The jet has settled into a quasi-steady state structure in the shape 
of a "corkscrew" (colors arranged from blue to red to represent 
low and high values of the density). The disk (not visible) is on 
the left-hand side of the image, and outflow is from left to right. 
[Adapted from OCP]. 

simulated jets maintain their long-term stability through a 
self-limiting process wherein the average Alfvenic Mach 
number within the jet is maintained to order unity. This 
is accomplished in at least two ways: (i) poloidal magnetic 
field is concentrated along the central axis of the jet form- 
ing a "backbone" in which the Alfven speed is sufficiently 
high to reduce the average jet Alfvenic Mach number to 
unity, and (ii) the onset of higher order Kelvin-Helmholtz 
"flute" modes (m > 2) reduce the efficiency with which 



the jet material is accelerated, and transfer kinetic energy of 
the outflow into the stretched, poloidal field lines of the dis- 
torted jet. This too has the effect of increasing the Alfven 
speed, and thus reducing the Alfvenic Mach number. The 
jet is able to survive the onset of the more destructive to = 1 
mode in this way. It was further discovered that jets go into 
alternating periods of low and high activity as the disappear- 
ance of unstable modes in the sub-Alfvenic regime enables 
another cycle of acceleration to super- Alfvenic speeds. The 
period of the episodic (low vs. high) behavior uncovered 
in the case of the wobbling jet (upper four panels in Fig.|3 
is of the order of 200 t^. For a typical young stellar object 
of mass 0.5Afo and for n ^ 0.05 AU (i.e. U - 0.9 day), 
this would correspond to a minimum period of roughly 180 
days. This is consistent with the temporal variations of the 
observed YSOs outflow velocities which appears to occur 
on timescales of a few years (Woitas et ai, 2002). 

Jets ultimately settle into relatively stable end-states in- 
volving either corkscrew (lower panel in Fig. |2j, or wob- 
bling types of structure (upper four panels in Fig. |2}. The 
difference in the two type of jets can be traced to the dif- 
ference in the profiles imposed at the accretion disk (see 
OCP for more details). This trend has been uncovered in 
3D simulations of large scale jets where wiggled jets have 
also been observed after the jets survived the destructive 
era (Uchida et al, 1992; Todo et al, 1993; Nakamura et 
al., 2001; Nakamura and Meier, 2004). For complete- 
ness, we should mention the recent similar simulations per- 
formed by Kigure and Shibata (2005) which included the 
disk self-consistently. These studies also show that non- 
axisymmetric perturbations dominate the dynamics. How- 
ever, these simulations were only run for very short times, 
never exceeding 2 inner-disk orbital periods due to numer- 
ical obstacles reported by the authors. In OCP simulations, 
it was shown that the Cartesian grid (also used by Kigure 
and Shibata) induces artificial modes in the early stages 
of the simulations and only after many orbital periods that 
these become washed out before the real modes enter the 
dynamics. Another suggested end-stated for 3D jets is a 
minimum-energy "Taylor state", wherein the jet maintains 
comparable poloidal and toroidal field components {Kdnigl 
and Choudhuri, 1985). This alternative does not pertain to 
our simplified simulations since we do not allow for internal 
magnetic energy dissipation. 

In summary, a large body of different simulations have 
converged to show that jets; (i) are centrifugally disk winds; 
(ii) are collimated by the "hoop stress" engendered by their 
toroidal fields; (iii) achieve two types of configuration de- 
pending on the mass loading that takes place in the un- 
derlying accretion disc - those that achieve collimation to- 
wards a cylinder, and those that have a wide-angle struc- 
ture; (iv) achieve two types of regimes depending on the 
mass loading that takes place in the underlying accretion 
disc - those that achieve a stationary state, and those that 
are episodic; (v) achieve stability by a combination of MHD 
mode coupling and "back-bone" effect leading to a self- 
regulating mechanism which saturates the instabilities; and 
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(vi) achieve different morphologies ranging from cylindri- 
cal wobbling to cork-screw structure, depending on the pro- 
file of imposed on the underlying accretion disk. 

4. COUPLED DISK- JET EVOLUTION 

While it is physically useful to regard the accretion disk 
as a boundary condition for the jet, this ignores critical is- 
sues such as the self-consistent radial distribution of the 
mass loading, or of the threading magnetic field across the 
disk. These, and other degrees of freedom can be found by 
solving the combined disk-outflow problem, to which we 
now turn. 

Significant progress on the launch of disk winds has been 
made by theoretical studies of a restricted class of self- 
similar 2D models (e.g. Wardle and Kdnigl, 1993; L/, 1995, 
1996; Ferreira, 1997; Casse and Ferreira, 2000). Since 
centrifugally driven outflows can occur for completely cold 
winds (BP82), the models have examined both "cold" and 
"warm" (i.e., with some kind of corona) conditions. In the 
former class, material from some height above the disk mid- 
plane must move upwards and eventually be accelerated 
outwards in a jet. Therefore, there must be a direct link 
between the physics of magnetized disks, and the origin of 
outflows. 

Calculations show that in order to match the properties 
of jets measured by Bacciotti et al. (2000), warm wind so- 
lutions are preferred, wherein a disk corona plays a central 
role. In this situation, gas pressure imbalance will assist 
with feeding the jet. A warm disk corona is expected on 
general grounds because it is a consequence of the magneto- 
rotational instability (MRI), as the vertically resolved disk 
simulations of Miller and Stone (2000) have shown. 

A semi-analytic model of the radial and vertical disk 
structure, that includes self-consistently the outflow, has 
also been presented by Campbell (2003). These solutions 
demonstrate explicitly that the outflow contributes to the 
loss of angular momentum in the disk by channeling it 
along field lines into the outflow. For self-consistent so- 
lutions to be possible, the turbulent Mach number has to be 
between 0.01 and 0.1. 

Models of outflows andjets generally assume ideal (non- 
resistive) MHD. However, inside the disk non-ideal effects 
must become important, because the accreting matter would 
otherwise never be able to get onto the field Unes that thread 
the disk and connect it with the outflow. In recent two- 
dimensional models, Casse and Keppens (2002, 2004) as- 
sumed a resistivity profile, analogously to the Shakura and 
Sunyaev prescription. This assumes the presence of an un- 
derlying turbulence within the disk. Only fairly large val- 
ues of the corresponding ass parameter of around 0.1 have 
been used. The subscript SS refers to Shakura and Sunyaev 
(1973), who were the first to introduce this viscosity param- 
eter. In all cases the system evolves to a steady equiUbrium. 
Using resistive simulations in a different context, Kuwabara 
et al. (2005) demonstrated that a substantial amount of 
energy can be transported by Poynting flux if the poloidal 



field falls off with distance no faster than r~^. Otherwise, 
the fast magnetosonic point is located closer to the Alfven 
point and the jet wiU be dominated by kinetic energy, which 
is the case in the simulations of Casse and Keppens (2004). 

The general stability of disk-outflow solutions is still 
being debated and the result may depend on the detailed 
assumptions about the model. In the solutions discussed 
here the accretion stress comes entirely from the large scale 
magnetic field rather than some small scale turbulence. As 
emphasized by Ferreira and Casse (2004), real disks have 
a turbulent viscosity just as they have turbulent magnetic 
diffusivity or resistivity. However, if the accretion stress 
does come entirely from the large scale magnetic field, the 
wind-driven accretion flows may be unstable {Lubow et al, 
1994). Kdnigl (2004) has shown recently that there are in 
fact two distinct solution branches-a stable and an unstable 
one. He argues that real disk/wind systems would corre- 
spond to the stable branch. 

Finally, the idea of a gently flared accretion disk is likely 
to be merely the result of the modeler's simplification rather 
than observational reahty. Indeed, accretion disks can be 
warped due to various instabilities which can be driven 
by radiation from the central object (e.g., Pringle, 1996), 
or, more hkely, by the outflow itself {Schandl and Meyer, 
1994). If a system is observed nearly edge-on, a warp in the 
accretion disk produces periodic modulations of the light 
curve. As Pinte and Menard (2004) have demonstrated, 
this may be the case in AA Tau. Observations frequently re- 
veal major asymmetries in bipolar outflows, which may be 
traced back to the corresponding asymmetries in the disk 
itself. The possible causes for these asymmetries may be 
either an externally imposed asymmetry such as one-sided 
heating by a nearby OB association, or an internal symme- 
try breaking of the disk- wind solution as a result increased 
rotation. Examples of the latter are familiar from the study 
of mean field dynamo solutions of accretion disks (Torkels- 
son and Brandenburg, 1994). 

4.1. Dead zones 

So far it has been assumed that the magnetic field is well 
coupled to the disk. This is certainly valid for most of the 
envelope of an accretion disk, which is ionized by a com- 
bination of cosmic rays, as well as X-rays from the central 
YSO. However, the deeper layers of the disk are strongly 
shielded from these ionizing agents and the degree of ion- 
ization plumets. 

This dense layer, which encompasses the bulk of the 
disk's column density, cannot maintain a sufficiently high 
electron fraction, and is referred to as the dead zone (Gam- 
mie, 1996). It is the poorly ionized region within which 
the MRI fails to grow as a consequence of the diffusivity of 
the field. Recent work of Fleming and Stone (2003) shows 
that, although the local Maxwell stress drops to negligible 
values in the dead zones, the Reynolds stress remains ap- 
proximately independent of height and never drops below 
approximately 10% of the maximum MaxweU stress, pro- 
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vided the column density in that zone is less that 10 times 
the column of the active layers. The non-dimensional ratio 
of stress to gas pressure is just the viscosity parameter, ass- 
Fleming and Stone (2003) find typical values of a few times 
10~^ in the dead zones and a few times 10~^ in the MRl- 
active layers. Inutsuka and Sana (2005) have questioned 
the very existence of dead zones, and proposed that the tur- 
bulent dissipation in the disk provides sufficient energy for 
the ionization. On the other hand, their calculation assumes 
a magnetic Reynolds number that is much smaller than ex- 
pected from the simulations (see Matsumura and Pudritz, 
2006; MP06). 

The radial extent of the dead zone depends primarily 
upon the disk column density as well as effects of grains and 
radiation chemistry (e.g., Sana et al, 2000; Matsumura and 
Pudritz, 2005, MP05; and MP06). Inside the dead zone the 
magnetic Reynolds number tends to be below a certain crit- 
ical value that is somewhere between 1 and 100 {Sana and 
Stone, 2002), making MRI-driven turbulence impossible. 
Estimates for the radial extent of the dead zone range from 
0.7-100 AU (Fromang et al., 2002), to 2-20 AU in calcu- 
lations by Semenov et al. (2004). For Chiang et al. (2001) 
models of disks that are well constrained by the observa- 
tions, and whose surface density declines as E oc Tq ', 
MP05 find that the extent of the dead zone is robust - ex- 
tending out to 15 AU typically, and is fairly independent of 
the ionizing environment of the disk. 

For smaller radii, thermal and UV ionization are mainly 
responsible for sustaining some degree of ionization. The 
significance of the reduced value of ass in the dead zones 
is that it provides a mechanism for stopping the inward mi- 
gration of Jupiter-sized planets (MP05, MP06). Jets can 
still be launched from the well-coupled surface layer above 
the dead zone (e.g., Li, 1996; Campbell, 2000). 

When the MRI is inactive in the body of the disk, alter- 
native mechanisms of angular momentum transport are still 
possible. In protostellar disks there are probably at least 
two other mechanisms that might contribute to the accre- 
tion torque: density waves (Rdzyczka and Spruit, 1993), and 
the interaction with other planets in the disk (Goodman and 
Rafikov, 2001). 

A more controversial alternative is to drive turbulence by 
a nonlinear finite ampUtude instability (Chagelishvili et al., 
2003). While Hawley et al. (1999) have given general ar- 
guments against this possibility, Afshordi et al. (2005) and 
other groups have continued invesligaling the so-called by- 
pass to turbulence. The basic idea is that successive strong 
transients can maintain a turbulent state in a continuously 
excited manner. Lesur and Longaretti (2005) have recently 
been able to quantify more precisely the critical Reynolds 
number required for instabihty. They have also highlighted 
the importance of pressure fluctuations that demonstrate 
that the general argument by Balbus et al. (1996) is insuffi- 
cient. 

Finally, the role of vertical (convectively stable) density 
stratification in disks and the possibiUty of the so-caUed 



strato-rotational instability has been proposed as a possi- 
ble mechanism for disk turbulence {Dubrulle et al, 2005). 
This instability was recently discovered by Molemaker et 
al. (2001) and the linear stability regime was analyzed by 
Shalybkov and Rudiger (2005). However, the presence of 
no-sUp radial boundary conditions that are relevant to ex- 
periments and used in simulations are vital. Indeed, the in- 
stability vanishes for an unbounded regime, making it irrel- 
evant for accretion disks ( Umurhan, 2006). 

4.2. Advected vs. dynamo generated magnetic fields 

One of the central unresolved issues in disk wind theory 
is the origin of the threading magnetic field. Is it dragged in 
by the gravitational collapse of an original magnetized core, 
or is it generated in situ by a disk-dynamo of some sort? 

The recent direct detection of a rather strong, true disk 
field of strength 1 kG at 0.05 AU in FU Ori, provides 
new and strong support for the disk wind mechanism (Do- 
nati et al. 2005). The observation technique uses high- 
efficiency high-resolution spectropolarimetry and holds 
much promise for further measurements. This work pro- 
vides excellent evidence that distinct fields exist in disks in 
spite of processes such as ambipolar diffusion which might 
be expected to reduce them. 

4.2.1 Turbulence effects. Accretion disks are turbulent 
in the well-coupled lower corona and beyond. They are 
also intrinsically three-dimensional and unsteady. Ques- 
tions regarding the stability of disks concern therefore only 
the mean (azimuthally averaged) state, where all turbulent 
eddies are averaged out. The averaged equations used to 
obtain such solutions incorporate a turbulent viscosity. In 
this cormection we must explain that the transition from a 
macroscopic viscosity to a turbulent one is more than just 
a change in coefficients, because one is really entering a 
new level of description that is uncertain in many respects. 
For example, in the related problem of magnetic diffusion it 
has been suggested that in MRI-driven turbulence, the func- 
tional form of of turbulent magnetic diffusion may be differ- 
ent, such that it operates mainly on the system scale and less 
efficiently on smaller scales (Brandenburg and Sokoloff, 
2002). It is therefore important to use direct simulations 
to investigate the stability in systems that might be unsta- 
ble according to a mean field description. There is no good 
example relevant to protostellar disks, but related experi- 
ence with radiatively dominated accretion disks, where it 
was found that the so-called Lightman-Eardly instability 
might not lead to the breakup of the disk (Turner et al, 
2004), should provide enough reason to treat the mean-field 
stability problem of protostellar disks with caution. Also, 
one cannot exclude that circulation patterns found in sim- 
ulations of the two-dimensional mean-field equations (e.g., 
Kley and Lin, 1992) may take a different form if the fully 
turbulent three-dimensional problem was solved. 

The turbulent regions of disks will generally be capable 
of dynamo action - wherein an ordered field is generated 
by feeding on the energy that sustains the turbulence. These 
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two aspects are in principle interlinked, as has been demon- 
strated some time ago using local shearing box simulations. 
These simulations show not only that the magnetic field is 
generated from the turbulence by dynamo action, but also 
that the turbulence itself is a consequence of the magnetic 
field via the MRI (Brandenburg et ai, 1995; Stone et al, 
1996). 

In comparing numerical dynamos with observational re- 
ahty, it is important to distinguish between two different 
types of dynamos: small scale and large scale dynamos (re- 
ferring primarily to the typical physical scale of the field). 
Both types of dynamos have in general a turbulent compo- 
nent, but large scale dynamos have an additional compo- 
nent on a scale that is larger than the typical scale of the 
turbulence. Physically, this can be caused by the effects of 
anisotropics, hehcity, and/or shear. These large scale dy- 
namos are amenable to mean field modeling (see below). 
Small scale dynamos, on the other hand, tend to be quite 
prominent in simulations, but they are not described by 
mean field models. However, small scale dynamos could be 
an artefact of using unreahstically large magnetic Prandtl 
numbers in the simulations. Protostellar disks have mag- 
netic Prandtl numbers, v / rj, around 10~®, so the viscos- 
ity of the gas, v, is much smaller than the diffusivity of 
the field, ry; see Table 1 of Brandenburg and Subrama- 
nian (2005). High magnetic Prandtl numbers imply that the 
field is advected with the flow without shpping too much, 
whereas for low numbers the field readily slips significantly 
with respect to the flow. At the resistive scale, where the 
small scale dynamo would operate fastest, the velocity is 
still in its inertial range where the spatial variation of the 
velocity is much rougher than for unit magnetic Prandtl 
number {Boldyrev and Cattaneo, 2004). This tends to in- 
hibit small scale dynamo action (Schekochihin etal., 2005). 
In many simulations, especially when a subgrid scale pre- 
scription is used, the magnetic Prandtl number is effectively 
close to unity. As a consequence, the production of small 
scale field may be exaggerated. It is therefore possible that 
in real disks large scale dynamo action is much more promi- 
nent than what is currently seen in simulations. 

In mean field models, only the large scale field is mod- 
eled. In addition to a turbulent magnetic diffusivity that 
emerges analogously to the turbulent viscosity mentioned 
above, there are also non-diffusive contributions such as the 
so-called a effect wherein the mean electromotive force can 
acquire a component parallel to the mean field, adynB; see 
Brandenburg and Subramanian (2005) for a recent review. 
For symmetry reasons, adyn has on average opposite signs 
above and below the midplane. Simulations of Ziegler and 
Riidiger (2000) confirm the unconventional negative sign of 
adyn in the upper disk plane, which was originally found in 
Brandenburg et al. (1995). This has implications for the 
expected field parity which is expected to be dipolar (anti- 
symmetric about the midplane); see Campbell et al. (1998) 
and Bardou et al. (2001). Although a tendency toward dipo- 
lar fields is now seen in some global accretion disk simula- 
tions {De Villiers et al., 2005), this issue cannot as yet be 



regarded as conclusive. 

4.2.2 Outflows from dynamo-active disks. Given that 
the disk is capable of dynamo action, what are the relative 
roles played by ambient and dynamo-generated fields? This 
question has so far only been studied in limiting cases. The 
first global simulations were axisymmetric, in which case 
dynamo action is impossible by the Cowling (1934) the- 
orem. Nevertheless, such simulations have demonstrated 
quite convincingly that a jet can be launched by winding 
up the ambient field and driving thereby a torsional Alfven 
wave (Matsumoto et al, 1996). These simulations are scale 
invariant and, although they were originally discussed in 
the context of active galaxies, after appropriate rescaling, 
the same physics applies also to stellar disk outflows. 

When the first three-dimensional simulations became 
available an immediate issue was the demonstration that 
MRI and dynamo action are still possible (Hawley, 2000). 
Although the simulations were applied to model black hole 
accretion disks, many aspects of such models are suffi- 
ciently generic that they carry over to protostellar disks as 
well. These simulations showed that the dynamo works effi- 
cientiy and produces normalized accretion stresses (ass ~ 
0.1) that exceed those from local simulations (ass ~ 0.01). 
The global simulations did not at first seem to show any 
signs of outflows, but this was mainly a matter of look- 
ing at suitable diagnostics, which emphasizes the tenuous 
outflows rather than the much denser disk dynamics (De 
Villiers et al., 2005). Although these simulations have no 
ambient field, a large scale field develops in the outflow re- 
gion away from the midplane. However, it is difficult to run 
such global simulations for long enough to exclude a de- 
pendence of the large scale field on the initial conditions. 
The outflows are found to be uncoUimated such that most 
of the mass flux occurs in a conical shell with half-opening 
angle of 25-30° (De Villiers et al , 2005). 

These very same properties are also shared by mean field 
simulations, where the three-dimensional turbulent dynam- 
ics is modeled using axisymmetric models {von Rekowski 
et al., 2003). In these simulations (Fig. 3), where the field 
is entirely dynamo-generated, there is an uncoUimated out- 
flow with most of the mass flux occurring in a conical 
shell with a half-opening angle of about 25-30°, just hke 
in the three-dimensional black hole simulations. The field 
strength generated by dynamoc action in the disk is found 
to scale as Bp oc r~^. This is very similar to the scahng for 
poloidal field strength found in collapse simulations (see 
Banerjee and Pudritz, 2006). 

The conical outflows discussed above have tentatively 
been compared with the observed conical outflows inferred 
for the BN/KL region in the Orion nebula out to a distance 
of 25-60 AU from its origin (Greenhill et al, 1998). One 
may wonder whether coUimated outflows are only possible 
when there is an ambient field (see Section 3). This conclu- 
sion might be consistent with observations of the Taurus- 
Auriga molecular cloud. Menard and Duchene (2004) 
found that, although T Tauri stars are oriented randomly 
with respect to the ambient field, there are no bright and ex- 
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Fig. 3. — Outflow from a dynamo active accretion disk 
driven by a combination of pressure driving and magneto- 
centrifugal acceleration. The extent of the domain is 
[0,8] X [—2,30] in nondimensional units (corresponding 
to about [0, 0.8] AU x [-0.2, 3] AU in dimensional units). 
Left panel: velocity vectors, poloidal magnetic field lines 
and gray scale representation of the specific enthalpy h in 
the inner part of the domain. Right panel: velocity vec- 
tors, poloidal magnetic field lines and normahzed specific 
enthalpy h/\^\ in the full domain. [Adapted from von 
Rekowski et al. (2003)]. 



tended outflows when the axis is perpendicular to the ambi- 
ent field. Note that Spruit and Uzdensky (2005) have argued 
that the efficiency of dragging an ambient field toward the 
star have been underestimated in the past, provided the field 
segregates into many isolated flux bundles. 

In summing up this section we note that all magnetic 
disk models produce outflows of some form. However, the 
degree of coUimation may depend on the possibility of an 
ambient field. An important ingredient that seems to have 
been ignored in all combined disk-jet models is the effect 
of radiation. There are also many similarities between the 
outflows in models with and without explicitly including the 
disk structure. Although these simulations are non-ideal (fi- 
nite viscosity and resistivity), the various Lagrangian invari- 
ants (mass loading parameter, angular velocity of magnetic 
field fines, as well as the angular momentum and Bernoulli 
constants) are stiU very nearly constant along field lines out- 
side the disk. 

5. JETS AND STAR-DISK INTERACTION 

In this section we review the recent development con- 
cerning the interaction of young stellar magnetospheres 
with the surrounding accretion disk and the contribution 



that these processes may have to the formation of jets. 

Perhaps the most important clue that star-disk interac- 
tion is an important physical process comes from the obser- 
vation that many YSOs are observed to spin at only a small 
fraction (~ 10%) of their break-up speed (e.g., Herbst et 
al, 2002). They must therefore undergo a significant spin- 
down torque that almost cancels the strong spin-up that aris- 
ing from gas accretion from the inner edge of the Keplerian 
disk. 

Three possible spin-down torques have been proposed: 
(i) disk-locking - a magnetospheric connection between the 
disk and the star which dumps accreted angular momen- 
tum back out into the disk beyond the co-rotation radius 
(Kdnigl 1991); (ii) an X-wind - a centrifugally driven out- 
flow launched from the very inner edge of the accretion 
disk which intercepts the angular momentum destined for 
the star (see chapter by Shang et al); and (iii) an accretion- 
powered stellar wind - in which the slowly rotating, cen- 
tral star drives a massive (Mw,* — O.lMa), stellar wind 
that carries off the bulk of the angular momentum as well 
as a fraction of the energy, that is deposited into the pho- 
tosphere by magnetospheric accretion {Matt and Pudritz, 
2005b). The last possibility implies that the central object 
in an accretion disk also drive an outflow, which is in ad- 
dition to the disk wind. Separating the disk wind from the 
outflow from the central object would be difficult, but the 
shear layer that separates them might generate instabihties 
and shocks that might be diagnostics. 

The strength of the star-disk coupling depends on the 
strength of the stellar magnetic field. Precise stellar mag- 
netic field measurements exist for nine TTSs, while for 
seven other stars statistically significant fields have been 
found (see Symington et al, 2005b and references therein). 
Zeeman circular-polarization measurements in He I of sev- 
eral classical TTs give direct evidence of kG magnetic fields 
with strong indication of considerable field variation on the 
time scales of years {Symington et al., 2005b). For a sam- 
ple of stars a magnetic field strength up to 4 kG fields could 
be derived (DF Tau, BP Tau), as for others, notably the jet 
source DG Tau, no significant field could be detected at the 
time of observation. 

5.1. Basic processes 

5.7.7 Variable accretion from a tipped dipole. A central 
dipolar field inclined to the rotation axis of star and disk 
may strongly disturbs the axisymmetry of the system. Pho- 
tometric and spectroscopic variability studies of AATau 
give evidence for time-dependent magnetospheric accretion 
on time scales of the order a month. Monte-Carlo modeling 
shows that the observed photo-polarimetric variability may 
arise by warping of the disk that is induced by a tipped mag- 
netic dipole {O'Sullivan et al, 2005). More general inves- 
tigations of the warping process by numerical simulations 
show that the warp could evolve into a steady state precess- 
ing rigidly {Pfeiffer and Lai, 2004). Disks can be warped by 
the magnetic torque that arises from the a slight misahgn- 
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Fig. 4. — Slice through a funnel stream. Density contours from 
p — 0.2 to p = 2.0. The density of the disk corona is p = 
0.01 — 0.02 Selected magnetic field lines. The rotational axis and 
the magnetic moment are indicated as fl and fi. [Adapted from 
Romanova et al. (2004).] 



ment between the disk and star's rotational axis {Lai 1999). 
Disk warping may also operate in the absence of a stellar 
magnetosphere since it can be induced by the interaction 
between a large-scale magnetic field that is anchored in the 
disk, and the disk electric current. This leads to a warp- 
ing instability and to the retrograde precession of magnetic 
jets/outflows {Lai, 2003). 

Three-dimensional radiative transfer models of the mag- 
netospheric emission line profile based on the warped disk 
density and velocity distribution obtained by numerical 
MHD simulations give gross agreement with observations 
with a variability somewhat larger than observed {Syming- 
ton et al, 2005a). 

5.1.2 Magnetic flux. Compared to the situation of a pure 
disk magnetic field, the magnetic field of the star may add 
substantial magnetic flux to the system. For a polar field 
strength Bq and a steUar radius i?*, the large scale stellar 
dipolar field (ignoring angular variations) 



Bp,^(r)~40G 



IkG/ VSE,^ 



(19) 



has to be compared to the accretion disk poloidal magnetic 
field which is provided either by a disk dynamo or by ad- 
vection of ambient intersteUar field. 
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where B^q is the equipartition field strength in the disk. 
This flux will not remain closed, but will inflate and open up 
as magnetic field lines are sheared and extract gravitational 
potential energy from the accreting flow (e.g., Uzdensky et 
al., 2002; Matt and Pudritz, 2005a). These field lines there- 
fore effectively become a disk field, and therefore follow 



the processes of disk wind production we have already dis- 
cussed. 

The additional Poynting flux that threads the disk may 
assist the jet launching by MHD forces and may serve as 
an additional energy reservoir for the kinetic energy of the 
jet implying greater asymptotic jet speed (Michel scaling; 
Michel, 1969; Fendtand Camenzind, 1996). 

5.7.5 Disk locking vs. stellar winds. The spin of the star 
will depend on how angular momentum arriving from the 
disk, is dealt with. In the magnetic "disk locking" picture, 
the threading field of the disk will re-arrange the global an- 
gular momentum budget. The torque on the star by the ac- 
cretion of disk matter is 



Ma (GM.nn) 
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(21) 



(e.g.. Matt and Pudritz, 2005a), with the disk accretion rate 
Maco the stellar mass and the disk inner radius Tin in- 
side the corotation radius Tco = (GMi)^/'^f7* On the 
other hand if "disk locking" is present, the stellar rotation 
may be decelerated by the magnetic torque due to stellar 
field lines connecting the star with the accretion disk out- 
side i?co- The differential magnetic torque acting on a disk 
annulus of dr width is 



rfTmag = r B^B^dr. 



(22) 



While Bz in principle follows from assuming a central dipo- 
lar field, the induction of toroidal magnetic fields is model 
dependent (disk resistivity, poloidal field structure). This is 
why recent numerical simulations of dipole-disk interaction 
that simultaneously evaluate the poloidal and toroidal field 
components have become extremely valuable. 

If indeed the star loses angular momentum to the disk 
(this is not yet decided by the simulations, see below), both 
disk accretion and jet formation are affected. In order to 
continue accretion, excess angular momentum has to be re- 
moved from the accreting matter. A disk jet can be an ef- 
ficient way to do this, as has been supposed in the X-wind 
picture. 

The central stellar magnetic field may launch a strong 
stellar wind to rid itself of the accreted angular momen- 
tum. Such an outflow will interact with the surrounding disk 
wind. If true, observed YSO jets may consist of two com- 
ponents - the stellar wind and the disk wind, with strength 
depending on intrinsic (yet unknown) parameters. The stel- 
lar wind (open field lines of stellar magnetosphere) exerts a 
spin-down torque upon the star of magnitude 



Mwind,*^*ri = 3 X 10' 
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As a historical remark we note that the model topology 
of dipole-plus-disk field were introduced for protosteUar jet 
formation more then 20 years ago in the MHD simulations 
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of Uchida and Shibata (1984). Further investigations con- 
sidering the detailed physical processes involved in disk 
truncation and channeling the matter along the dipolar field 
Unes by Camenzind (1990), Kdnigl (1991) and Shu et al. 
(1994) resulted in a breakthrough of these ideas to the pro- 
tostellar jet community. 

5.2. Numerical simulations of star-disk interaction 

The numerical simulation of the magnetospheric star- 
disk interaction is technically most demanding since one 
must treat a complex model geometry in combination with 
strong gradients in magnetic field strength, density and re- 
sistivity. In general, this may imply a large variation in 
physical times scales for the three components of disk, jet, 
and magnetosphere, which all have to be resolved numeri- 
cally. Essentially, numerical modeling of the star-disk inter- 
action requires a diffusive and viscous MHD code including 
radiative transfer. In addition, realistic models for recon- 
nection processes and disk opacity are needed. Then, such 
a code has to run with high spatial and temporal resolution 
on a global scale. 

Compared to the situation about a decade ago when there 
was still only Uchida and Shibata's (1984) initial (though 
ten years old) simulation available, today huge progress has 
been made with several groups (and also codes) competing 
in the field. Early simulations were able to follow the evolu- 
tion only for a few rotations of the iimer disk (note that 100 
rotations at the co-rotation radius correspond to 0.3 Keple- 
rian rotations at 10 co-rotation radii) (Hayashi et al., 1996; 
Miller and Stone, 1997; Goodson et al, 1997). Among the 
problems involved is the initial condition of the simulation, 
in particular the nature of the disk model which could be 
numerically treated. A steady initial corona will strongly 
shear with the Keplerian disk leading to current sheet (thus 
pressure gradients) along the disk surface. If the simula- 
tions run long enough, this could be an intermittent feature. 
However, the danger exists that the artificial current sheet 
will fatally destroy the result of the simulation. Applying 
a Shakura-Sunyaev disk for the initial disk structure, the 
code should also consider 0:55 viscosity. Otherwise the ini- 
tial disk evolution is not self-consistent (see, e.g., Goodson 
etal, 1997). 

The next step is to increase the grid resolution and to 
redo the axisymmetric simulations. Goodson et al. (1999) 
and Goodson and Winglee (1999) were able to tteat sev- 
eral hundreds of (inner) disk rotations on a global grid of 
50 AU extension. The main result of these simulations is a 
two component flow consisting of a fast and narrow axial 
jet and a slow disk wind, both launched in the inner part 
of the disk (;• < 1 AU). An interesting feature is that the 
narrow axial jet is actually a coUimation in density and not 
in velocity. Close to the inner disk radius repetitive recon- 
nection processes are seen on time scales of a couple of 
rotation periods. The dipolar field inflates and an expand- 
ing current sheet builds up. After field recormection along 
the current sheet, the process starts aU over again. The os- 



cillatory behaviour leads to the ejection of axial knots. On 
average the central dipolar magnetosphere remains present 
and loaded with plasma. Forbidden emission fine intensity 
maps of these simulations have been calculated and allow 
for direct comparison with the observations. However, the 
numericaUy derived time and length scales of axial knots 
were still too different from what is observed. The magne- 
tospheric origin of jets (stellar dynamo) in favor of a pure 
disk origin ("primordial field") has also been stressed by 
Matt etal. (2002). 

In order to be able to perform long-term simulations of 
dipolar magnetospheres interacting with the accretion disk, 
Fendt and Elstner (1999, 2000) neglected the evolution of 
the disk structure and instead assume that the disk is a 
fixed boundary condition for the outflow (as in OPI). Af- 
ter 2000 rotations a quasi-steady state was obtained with 
a two-component outflow from disk and star. The outflow 
expands almost radially without signature of colhmation on 
the spatial scale investigated (20 x 20 inner disk radii). One 
consequence of this very long simulation is that the axial 
narrow jet observed in other simulations might be an inter- 
mittent feature launched in the early phase of the simulation 
as the initial field is reconfigured to a new equilibrium. The 
axial outflow in this simulations is massive but slow, but 
tends to develop axial instabilities for a lower mass loading. 
Clearly, since the disk structure was not been taken into ac- 
count, nothing could be said about the launching process of 
the outflow out of the disk. 

In a series of ideal MHD simulations Romanova and col- 
laborators succeeded in working out a detailed and suffi- 
ciently stable numerical model of magnetospheric disk in- 
teraction. They were the first to simulate the axisymmetric 
funnel flow from the disk inner radius onto the stellar sur- 
face (Romanova et al., 2002) on a global scale (-Rmax = 
50i?in) and for a sufficiently long period of time (300 rota- 
tions) in order to reach a steady state in the accretion fun- 
nel. The funnel flow with free-falling material builds up as a 
gap opens up between disk and star. The authors further in- 
vestigated the angular momentum exchange due to the disk 
"locking" by the funnel flow. Slowly rotating stars seem to 
break the inner disk region and spin up, while rapid rotators 
would accelerate the inner disk region to super-Keplerian 
velocity and slow down themselves. Only for certain stel- 
lar rotational periods a "torque-less" accretion could be ob- 
served. Strong outflows have not been observed for the 
parameter space investigated, probably due to the matter 
dominated corona which does not allow for opening-up the 
dipolar field. 

Further progress has been achieved extending these sim- 
ulations to three dimensions (Romanova et al, 2003, 2004). 
For the first time it has been possible to investigate the inter- 
action of an inclined stellar dipolar magnetosphere with the 
surrounding disk. For zero inclination axisymmetric fun- 
nels are found as in the axisymmetric simulations. For non- 
zero inclination the accretion splits in two main streams fol- 
lowing the closest path from the disk to the stellar surface as 
is seen in Fig. 4. Magnetic braking changes the disk struc- 
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ture in several ways - its density structure (gaps, rings), its 
velocity structure (sub-Keplerian region) and its geometry 
(warping). The slowly rotating star is spun up by accret- 
ing matter but this acceleration does only weakly depend 
on inchnation. For completeness we should note that for 
these simulations a clever approach for the numerical grid 
has been applied, the "cubed grid" which does not obey the 
singular axes as for spherical coordinates. 

The star-disk couphng by the stellar magnetosphere was 
also investigated by Kiiker et al. (2003). These simulations 
have been performed in axisymmetry, but an advanced disk 
model has been applied. Taking into account a- viscosity, a 
corresponding eddy magnetic diffusivity and radiative en- 
ergy transport, the code enables the authors to treat a reahs- 
tic accretion disk in their MHD simulations. As a result the 
simulations could be advanced to very long physical time 
steps into a quasi stationary state of the disk evolution. The 
authors show that the commonly assumed 1000 G magne- 
tosphere is not sufficient to open up a gap in the disk. Un- 
steady outflows may be launched outside of the corotation 
radius with mass loss rates of about 10 % of the accretion 
rate. The authors note that they were unable to detect the 
narrow axial jet seen in other publications before (in agree- 
ment with Fendt and Elstner, 2000). Magnetic braking of 
the star may happen in some cases, but is overwhelmed 
by the accretion torque still spinning up the star (however, 
these simulations treat the stellar surface as inner boundary 
condition). 

Using a different approach, the boundary condition of 
a stellar dipolar magnetosphere or a large scale disk field 
were relaxed. This allows for the interaction of a stellar 
dipole with a dynamo-generated disk magnetic field (von 
Rekowski and Brandenburg, 2004), or the evolution of a 
dynamo-generated stellar field affected by a surrounding 
disk with its own disk dynamo (von Rekowski and Bran- 
denburg, 2006) to be studied. The results of this work 
are in agreement with previous studies, and directly prove 
that a disk surrounding a stellar magnetosphere may actu- 
ally develop its own magnetic field strong enough to launch 
MHD outflows. In agreement with results by Goodson and 
Winglee (1999), accretion tends to be unsteady and alter- 
nating between connected and disconnected states. Even 
for reahstic stellar fields of several hundred gauss the mag- 
netic spin-down torque is insufficient to overcome the spin- 
up torque from the accretion flow itself. As shown by 
the authors, angular momentum exchange is complex and 
may vary in sign along the stellar surface, braking some 
parts of the star and accelerating others. The dynamo- 
generated stellar field may reach up to 750 G and switches 
in time between dipolar and quadrupolar symmetry, while 
the dynamo-generated 50 G disk field is of dipolar symme- 
try. In general, the dynamo-generated stellar field is better 
suited to drive a stellar wind. The simulations show that for 
these cases stellar wind braking is more efficient than brak- 
ing by the star-disk coupling, confirming the results of Matt 
and Pudritz (2005a). 

Recent studies of torque-less accretion (Long et al. 



2005) compare cases of a (i) weak stellar magnetic field 
within a dense corona and (ii) a strong field in a lower den- 
sity corona. They investigate the role of quasi-periodic field 
line reconnection in couphng the disk and stellar fields. Un- 
like previous works, these authors conclude that magnetic 
interaction effectively locks the stellar rotation to about 
10% the break-up velocity, a value which actually depends 
on the disk accretion rate and the stellar magnetic moment. 
While they correctly stress the importance of deahng with 
the exact balance between open and closed field Unes and 
their corresponding angular momentum flux, the exact mag- 
netospheric structure in the innermost region will certainly 
depend on the resistivity (diffusivity) of the disk and corona 
material. This is, however, not included in their treatment 
as an ideal MHD code has been apphed. 

We summarize this section noting that tremendous 
progress has occurred in the numerical simulation of the 
star-disk interaction. The role of numerical simulations 
is pivotal in this area because the mechanisms involved 
are complex, strongly interrelated, and often highly time- 
dependent. It is fair to say that numerical simulations of the 
star-disk interaction have not yet shown the launching of a 
jet flow comparable to the observations. 

6. JETS AND GRAVITATIONAL COLLAPSE 

Jets are expected to be associated with gravitational col- 
lapse because disks are the result of the collapse of rotating 
molecular cloud cores. One of the first simulations to show 
how jets arise during gravitational collapse is the work of 
Tomisaka (1998, 2002). Here, the collapse of a magnetized 
core within a rotating cylinder of gas gave rise the formation 
of a disk from which a centrifugally drive, disk wind was 
produced. Analytical work by Krasnopolsky and Konigl 
(2002) provides an interesting class of self-similar solutions 
for the collapse of rotating, magnetized, isothermal cores 
including ambipolar diffusion that could be matched onto 
BP82 outflows. 

Given the importance of Bonner-Ebert (B-E) spheres in 
the physics of star formation and gravitational collapse, re- 
cent efforts have focused on understanding the evolution of 
magnetized B-E spheres. Whereas purely hydrodynamic 
collapses of such objects never show outflows (e.g., Fos- 
ter and Chevalier, 1993; Banerjee et al, 2004), the addi- 
tion of a magnetic field produces them. Matsumoto and 
Tomisaka (2004) have studied the collapse of rotating B-E 
spheres wherein the magnetic axis is inclined with the ini- 
tial rotation axis. They observe that after the formation of an 
adiabatic core, outflow is ejected along the local magnetic 
field lines. Eventually, strong torques result in the rapid 
aUgnment of the rotation and magnetic field vectors. 

The collapse of a magnetized, rotating, B-E sphere with 
molecular cooling included (but not dust), was carried out 
by Banerjee and Pudritz (2006) using the FLASH adaptive 
mesh refinement code. The results of this simulation are 
shown in Fig. 5 which shows the end state (at about 70,000 
yrs) of the collapse of a B-E sphere that is chosen to be 
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Fig. 5. — Large scale outflow (top frame, scale of hundreds 
of AU), and and small scale disk wind and jet formed (bot- 
tom frame, scale of a fraction of an AU) during the gravi- 
tational collapse of a magnetized B-E, rotating cloud core. 
Cross-sections through the disk and outflows are shown - 
the blue contour marks the Alfven surface. Snapshots taken 
of an Adaptive Mesh calculation at about 70,000 years into 
the collapse. [Adapted from Banerjee and Pudritz, 2006]. 

precisely the Bok globule observed by Alves et al. (2001) - 
whose mass is 2.1Mq and radius R = 1.25 x lO'^ AU at 
an initial temperature of 16 K. Two types of outflow can be 
seen: (i) an outflow that originates at scale of ~ 130 AU 
on the forming disk that consists of a wound up column of 
toroidal magnetic field whose pressure gradient pushes out 
a slow outflow; and (ii) a disk-wind that collimates into a 
jet on scale of 0.07 AU. A tight proto-binary system has 
formed in this simulation, whose masses are still very small 
< 1O~^M0, which is much less than the mass of the disk 
at this time ~ IO'^Mq. The outer flow bears the hallmark 



of a magnetic tower, first observed by Uchida and Shibata, 
and studied by Lynden-Bell (2003). Both flow components 
are unbound, with the disk wind reaching 3 km s^^ at 0.4 
AU which is quite super-Alfvenic and above the local es- 
cape speed. The outflow and jet speeds will increase as the 
central mass grows. 

We conclude that the theory and computation of jets and 
outflows is in excellent agreement with new observations 
of many kinds. Disk winds are triggered during magne- 
tized collapse and persist throughout the evolution of the 
disk. They efficiently tap accretion power, transport a sig- 
nificant part portion of the disk's angular momentum, and 
can achieve different degrees of collimation depending on 
their mass loading. Accretion-powered stellar winds may 
also solve the stellar angular momentum problem. We are 
optimistic that observations will soon be able to test the 
universality of outflows, all the way from circumplanetary 
disks to those around O stars. 
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